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Abstract A review of various exactly solvable models on the determination of 
the energy spectra E{k) of 3D-velocity field, induced by chaotic vortex lines is 
proposed. This problem is closely related to the sacramental question whether a 
chaotic set of vortex filaments can mimic the real hydrodynamic turbulence. The 
quantity < v(li)v(— k) > can be exactly calculated, provided that we know the 
probability distribution functional 3^{{s{^,t)}) of vortex loops configurations. 
The knowledge of ^({s(^,f)}) is identical to the full solution of the problem of 
quantum turbulence and, in general, ^ is unknown. In the paper we discuss sev- 
eral models allowing to evaluate spectra in the explicit form. This cases include 
standard vortex configurations such as a straight line, vortex array and ring. In- 
dependent chaotic loops of various fractal dimension as well as interacting loops 
in the thermodynamic equilibrium also permit an analytical solution. We also de- 
scribe the method of an obtaining the 3D velocity spectrum induced by the straight 
line perturbed with chaotic ID Kelvin waves on it. 



PACS number(s): 67.25.dk, 47.37.-i-q, 03.75.Kk 



1 Introduction 

One of the exciting applications of quantum turbulence is the solution (or rather, 
the attempt at a solution) of the tantalizing problem of classical turbulence. For 
this reason, the problem of the spectrum of the 3D velocity field induced by a 
chaotic vortex filament becomes one of the central questions. The formal rela- 
tion, allowing of calculating £(k) = Ps (vkV_k) via the vortex line configuration 
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{s(^)}, is given by formula (se©^) 

/ -:. ^ L L 



2 J /s'(^i)s'(^2)^^i^^2exp 
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In isotropic case, the spectral density depends on the absolute value of the wave 
number k. Integrating over solid angle leads to formula (seOi^): 



i7(i.\ I P'^^ ) f.'d^ \.'(^ sin(^|s(^i)-s(fe)|) \ 
= \ (2^/ [ ' ^^^^^^^^^^ ^|s(4i)-s(&)| / • 

For anisotropic situations, formula (|2]i is understood as an angle average. Further 
we will apply these formulas to study some particular situation. 
Of course, the relations ©-([Tl are just mathematical identities, and the physics is 
hidden behind the () operation. To calculate this average, we need the probability 
distribution functional ^({s(^)},f), which is the probability that the system has 
the configuration {s(^)} of the set of the vortex loops. Knowing ,^({s(^)},f) 
is identical to the full solution of the problem of quantum turbulence. Therefore, 
the problem of the determination of the 3D energy spectrum is not resolved, al- 
though there are a series of theoretical approaches and numerical simulations. We 
consider here several exactly solvable cases, the study of them is quite instructive. 



2 Regular structures (straight line, vortex array, ring) 

The angle averaged spectrum created by straight vortex line (directed along axis 
Z) is directly evaluated from (|2| 

k \0 ' k 

The length L is the conditional quantity, more correctly to talk here about energy 
per unit length. Spectrum of the straight line ([3]) is discussed early (see, e.g., »,—,), 
Vinen^ proposed the spectrum on the basis of dimensional consideration. This 
result is important, since it ([3]) states that for any vortex system the high wave 
numbers larger than inverse curvature, E(]d) should be (ps"^ l\n)k^^ (per unit 
length). 

Let us take a set of straight vortex filaments forming the square lattice U ) = 
U('t(,)'/,z)- Points x,,y, are coordinates for vortices on the xy-plane, index / runs 
from 1 to A'^ . The neighboring lines are separated by distance h. I.e., — O. 

Then the general relation OJ leads to the following formula (Hanninen (private 
communication), Nowak et al. (^): 
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where dij = ^ (x/ — Jf/)^ + iji — J/)^ distances between vortices on the xy-plane. 
Thus, determination of the spectrum on the basis (|4]l should be done with the use 
of the quadruple summation (over (x,,Xy,y,,yy)), which requires large computing 
resources. Clear, however, that for very small k, which corresponds to very large 
distance, the whole array can be considered as large single vortex with the circula- 
tion A^^k:. Accordingly, the spectrum (per unit height) should be (pjA'^K:^/47r)fc~'. 
For large k, which corresponds to very small distance from each line, the spectrum 
(per unit height) should be (PsK:^/47r)fc^' In the intermediate region kb « 1, and 
Nkb >> 1 (this condition implies that inverse wave number k~^ is larger intervor- 
tex space between neighboring lines, but smaller then the size of the whole array 
Ndij), we can replace the quadruple summation by the quadruple integration with 
infinite limits. This procedure corresponds that we exclude the fine-scale motion 
near each of vortex, and are interested in the only large-scale, coarse-grained mo- 
tion. After obvious change of variables — ?> fcx/, — ?> ky, etc. we get that the the 
whole integral should scale as 1 /k'^, and accordingly 

E{k)o:\/k^ (5) 

Formula E (k) = dE jdk oc 1 jk^ implies that Ei^=dE jSk should behave as 
1 lk~^. It follows, for example, from that in the 2D case d^k = Inkdk (isotropic 
is assumed). Since, further, ii(k) = v(k)v(— k), we see that v(k) scales as 1 /k~^ . 
The latter means (see, e.g.,^, Eqs. (4.60), (4.61)) that the velocity v(r) scales as r'. 
Thus, the uniform vortex array creates the course-grained motion, which is rota- 
tion, as it should be. Moreover, the coefficient is equal to K/2b^, which coincides 
with the Feynman rule. 

Now we will consider a vortex ring with radius R lying in the x-y plane. The 
line s{^) can be parameterized as s(x) = (7?cos(p,7?sin(p,0) with (p G [0,27r]. 
Applying it to Eq. (|2]l we get (see also^) 



E, 



^K^R f f COS 91-92) sin(2^j?sin (pi-(i(?2)/2 ) 

-^JJdcp^dcp, 2fcsm((<p.-9z)/2) 



Evaluating the integral numerically shows that the spectrum Eri„g{k) scales 
like Erini, [k) ~ k^ for kR<^\. Frequently, the k^ is referred to as a proof for the 
thermodynamical equilibrium state. We would like to stress, however, that this 
distribution of the energy (valid far from the ring) has nothing to do with the 
equipartition law. It is a consequence of the fact that closed vortex loops induce a 
far field flow scaling as l/P . That, in turn, generates a spectrum E{k) oc k^. This 
fact was established for classical turbulence (see, e.g.,**). For quantum turbulence 
this result was discussed by Stalp, Skrbek and Donnelly—. For the large k, namely 
for ^ 1, spectrum E,if,^{k) scales like ^ k~^ as for straight line. 



ring - fr,.2 



3 Noninteracting Gaussian loops of various fractal dimension 

One of the approaches, allowing exact solution, based on the viewing of the VT 
as a set of loops having the random structure with the various fractal dimension ^"i. 
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Fig. 1 (Color online) Energy spectra of a pure fractal line (left), and a more realistic semi-fractal 
vortex loop (right). 



This theory is based on the Gaussian model (seeii). The probability ^Gmiss{{s{^ : }) 
of finding a particular configuration {s(^ ,f)} is expressed by the probability dis- 
tribution functional (for details, see the paper by author—) 

(7) 

Here ^/K is a normalizing factor and I is the length of curve. The typical form of 
function — is a smoothed 5 function, having a Mexican hat shape with 

width equal to ^o- Behavior of A^jii^i — ^2) on large scales determines the sec- 
ond order correlation function between tangent vectors (s'„(^i)s'„(^2)); which, in 
turn, determines the fractal dimension of line. Referring the reader to the original 
article for the details we briefly state the main results. 

In the region ^ ^ 1 /i'/^^, the spectrum behaves asE{k) k^. We again stress 
that this dependence is related to the far field flow 1 /r^ from the restricted domain 
of vorticity (but not to the thermodynamical equilibrium). Note, that L}!^'^ is noth- 
ing but the real 3D size of the loop. In the opposite case, of large k we get 

This result had been obtained earlier from qualitative considerations and is dis- 
cussed in Frisch^. Fig. [T| (left) depicts (in logarithmic scale) three curves, the 
spectral densities E{k) for a pure fractal vortex filament of length L = 100 (all 
units are arbitrary) and three different values of Hd = 1,5/3,2 correspondingly. 
It can be seen how E{k) 0= is changed to a dependence E{k) 0= k^'^^^'^ in the 
region of large ^ 1/100" . 

Vortex loops in superfluid turbulent quantum fluids are "semifractal" objects, 
which are smooth at small — £,") along the curve, and are fully uncorrelated 
for the remote parts. A good approximation for (Sa(^i)Sa(^2)) is a function of the 
type 1/(1 + (4^/^0)2). 

Fig- [H (right) shows (in logarithmic scale) the spectral densities E{k) for a 
"semifractal" vortex loop of length L = 100 with^o=0.01 and(^o= 1-Itis clearly 
seen that the curves (^0 = 0.01) have three regions with bends at the points ki = 
1 / \/LE,q = 1, and kr = 1 / ^0 = 100. In these three different regions we have E{k) o< 
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k^, E [k) oc k^ and E{k)^k \ correspondingly. This is in excellent agreement with 
the qualitative considerations. 



4 Thermodynamically equilibrium vortex tangle 

In the previous Section, we considered the model of free (norunteracting) vortex 
loops. Let us now discuss another solvable case, that of interacting loops in ther- 
modynamical equilibrium. This can be accomplished by introducing the Boltz- 
mann factor exp(— j3// {s})into the Gaussian model O (see ^'■^^). Then the parti- 
tion function can be written (we will take here the purely Wiener distribution with 
A being the 5 function and with an elementary step equal to ^o) 







Z= / Ds(^}exp -A /(s'(^))2^^ -/^^ / i d^^d^2-^^^^^^,d^id^2 
J 2^qJ %n J J |s(^i) -s(^2)r 

(9) 

Here j3 = 1 /ksT is the inverse temperature. A considerable simplification in 
the evaluation of the partition function can be reached with the use of the Edwards 
trick (see for details^,— ,—,—), namely, exp(— j3// {s}) can be written as a the 
Gaussian path integral over an auxiliary vector field A(r). After all the transfor- 
mations, the partition function acquires the form of a Gaussian path integral over 
the 3D auxiliary vector field A(r) and E = dlnZ/dfi is (see ^^): 

E = ^^oLp...7..^^ (10) 

with = P^oLpsK^ /2y (Y is the volume). For small k, the energy spectrum 
E{k) cc {kBT)k^, which is nothing but the Rayleigh equipartition law dE/d^k =const. 
For large wave numbers E{k) is just about constant. The reason is that for large k, 
the interaction energy in the partition function (|9]) is smaller and the main contri- 
bution appears from the configuration term (s'(^))^, which is related to the con- 
nectivity of the line. Hence, we have the same spectrum as for a pure random walk 
with Hd = 2, discussed in previous Section. 

5 ID Kelvin waves spectrum and 3D velocity spectrum 

In the literature there is discussed the idea of obtaining the 3D velocity spectrum 
just by putting it equal to the spectrum of ID Kelvin waves. For instance, as stated 
inii:' 'We notice that, because the fluctuations of the velocity field are induced by 
the Kelvin wave fluctuations on the filaments, it is reasonable to expect that 

E{k)^EKw{ky\ (11) 

The same idea was used in papers by L'vov et al.(see e.g., --). Details of this activ- 
ity can be read in a series of papers by L'vov, Nazarenko and coauthors i2.iSi22i^ 
Let us consider this problem on the basis of general formula ([Tl (see'^). We 
takes(^ ,f) = {x{z^t)^y{z,t),z) and denote the two-dimensional vector {x{z^t),y{z,t)) 
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as ap (z, f ) (where the dimensionless amplitude a ^ 1). Substituting it into (|2]| and 
expanding in powers of a, we get, 

sin(fc|z2-zi|)(p(z2)-p(zi))^ ^ {p'{zi)p'{z2)sin{k\z2-zi\) ^^^^ 
2k\z2-zi\^ k\z2-zi\ 

The first term Eq of the zero-order in amplitude a exactly coincides with the energy 
spectrum induced by the unperturbed straight line (O, as it should be. 

To move further we have to find the correlation characteristics for the fluctu- 
ating vector of displacement p{z2)- We accept that the ensemble of Kelvin waves 
has a following power-like spectrum: 



{p{p)p{-p))=Ap-\ (14) 

We take here the notation p for the one-dimensional vector, conjugated to z, re- 
serving the notation k for the absolute value of the wave vector of the 3D field. The 
formula (114b implies that (see, e.g.,^, Eqs. (4.60), (4.61)) the squared increment for 
the vector of displacement scales as, ((p(z2) ^P(zi))^) cc {z2 — Then the 
second order correlator ((p'(z2)p'(zi))) scales as ((p'(z2)p'(zi))) cc {z2~ziY~^. 
Substituting it into (1141) and counting the powers of quantity k, we conclude that 
the correction dE{k) to the spectrum E{k), due to the ensemble of Kelvin waves 
has a form: 

dE{k)o:a^k-'+^. (15) 

It is remarkable fact that this quantity coincides formally with the one-dimensional 
spectrum of KW 5E{p) cc a^p~^'^^ however this contribution is small, by virtue 
the smallness of the wave amplitudes a, and disappears with the KW. 



6 Conclusions 

Summarizing, it can be concluded that the 3D energy spectrum E{k) consists of 
several parts. At small k, associated with the large scales, on the order of the size 
of the system, or on the scale of the stirring forcing (grids, propellers, vibrating 
objects, etc.), the spectrum behaves as E{k) oc k^. We recall again that this is 
the consequence of the asymptotic behavior {r — > oo) of the velocity field, not 
of thermodynamical equilibrium, for large wave numbers, exceeding the inverse 
intervortex space k > 2n/5, the energy spectrum E{k) should be close to k^^, 
again regardless of the specific model. The region of intermediate k is the most 
intriguing and exciting, it depends on fractal properties of lines. Unfortunately, so 
far there is no theory, predicting the Kolmogorov like spectrum ii(fc) cc k^^^^. 
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